First-Matsubara-frequency rule in a Fermi liquid. 
Part II: Optical conductivity and comparison to experiment 
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Motivated by recent optical measurements on a number of strongly correlated electron systems, 
we revisit the dependence of the conductivity of a Fermi liquid, <r(f2,T), on the frequency fi and 
temperature T. Using the Kubo formalism and taking full account of vertex corrections, we show 
that the Fermi liquid form Re<r -1 (f2,T) oc fl 2 + 4n 2 T 2 holds under very general conditions, namely 
in any dimensionality above one, for a Fermi surface of an arbitrary shape (but away from nesting 
and van Hove singularities), and to any order in the electron-electron interaction. We also show that 
the scaling form of Recr _1 (Sl, T) is determined by the analytic properties of the conductivity along 
the Matsubara axis. If a system contains not only itinerant electrons but also localized degrees of 
freedom which scatter electrons elastically, e.g., magnetic moments or resonant levels, the scaling 
form changes to Recr _1 (fi, T) oc Q. 2 + biv 2 T 2 , with 1 < b < oo. For purely elastic scattering, 6=1. 
Our analysis implies that the value of 6 ~ 1, reported for URu2Si2 and some rare-earth based doped 
Mott insulators, indicates that the optical conductivity in these materials is controlled by an elastic 
scattering mechanism, whereas the values of 6 ~ 2.3 and 6 ~ 5.6, reported for underdoped cuprates 
and organics, correspondingly, imply that both elastic and inelastic mechanisms contribute to the 
optical conductivity. 

PACS numbers: 71.10.Ay, 71.10. Pm 



I. INTRODUCTION 

Optical response of strongly correlated materials is an 
invaluable tool for studying the dynamics of charge car- 
riers i On par with the angular-resolved photoemission 
spectroscopy which gives information about the single- 
particle self-energy, optical experiments provide informa- 
tion about two other important quantities: the dynam- 
ical effective mass and the scattering rate of conduction 
electrons. 

In the preceding paper 2 (hereafter referred to as I), we 
discussed constraints imposed on the functional form of 
the retarded single-particle self-energy T, R (co,T) by the 
"first-Matsubara-frequency rule" . This rule stipulates 
that, under certain conditions, a function obtained by an- 
alytic continuation of Im£ (a>,T) to the Matsubara axis 
must vanish at the first fermionic Matsubara frequency 
uj — > ztinT. The familiar scaling form of ImE fl (w, T) in 
a generic Fermi liquid (FL) in D > 2 



ImE R {u,T) = C(u 



7T 2 T 2 ) 



(1.1) 



(with a coefficient of exactly 7r 2 in front of the T 2 term), 
obviously satisfies this rule. 

In the present paper, we discuss a similar constraint- 
"the first bosonic Matsubara-frequency rule"- imposed 
on the scaling form of the optical conductivity er(il,T). 

Within the semiclassical Boltzmann equation, the 
T 2 -scaling of the dc resistivity due to Umklapp 
electron-electron scattering was obtained by Landau and 
Pomeranchuk, 3 and due normal scattering in a two-band 
metal by Baberji Later on, Eliashberg^ re-derived this 
result from the Kubo formula, and showed that it re- 
mains valid to all orders in the electron-electron interac- 
tion. The ft/T scaling of the "optical resistivity" of a 3D 



FL was first discussed by Gurzhi^ who used a quantum 
Boltzmann equation to show that 

Re[p(0,r)] = Re[CT~ 1 (ri,T)] = A' [Q 2 +4ir 2 T 2 ] . 

(1.2) 

The f2 and T dependences of Re p(Q,, T) are similar to 
those of the leading term in ImE R (w,T) [cf. Eq. (fTTT]) ]. 
but the ratio of the T 2 and ft 2 terms is now 4-7T 2 instead 
of tt 2 . 

This difference is not accidental. Indeed, ImT, R (ui,T) 
measures the decay rate of single-particle excita- 
tions, which are fermions; hence the thermal part of 
ImE fl (oj,T) contains the square of first fermionic Mat- 
subara frequency (= -kT) rather than T itself. On the 
other hand, Rep(51, T) measures the decay rate of cur- 
rent fluctuations, which are bosons; hence the thermal 
part of Recr(f2,T) contains the square of the first (non- 
zero) bosonic Matsubara frequency (= 2ttT). Also not 
coincidentally, Eq. (|1.2p is of the same form as the sound 
absorption rate in FLs. 7 

To the best of our knowledge, the scaling form pre- 
dicted by Eq. (|1.2p has never been verified experimen- 
tally in conventional metals. On the other hand, the 
fi/T scaling of Rep(51, T) has been studied intensively in 
strongly correlated materials, e.g., in heavy- fermion met- 
als and doped Mott insulators. The result of these studies 
is quite surprising: whenever it was possible to fit the f2 
and T dependencies of Rep(f2, T) by quadratic functions, 
the coefficient b = 7T 2 T 2 /f2 2 was found to be quite differ- 
ent from 4. This issue was highlighted by recent study 8 
of Rep(f2, T) in the "hidden-order" (HO) heavy-fermion 
compound URii2Si2, where b was found to be close to 1 
above the 17.5 K transition to the HO state. In fact, the 
value of b w 1 was found in a number of other materials, 
including two rare-earth based doped Mott insulators, 
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NcWTiOa (Ref.i) and Ceo.095Cao.05TiO3.04 (Ref.HS). 
Another recent study-ii reports b « 2.3 in the underdoped 
cuprate HgBa 2 Cu04 + 5. Whereas the observed value of 
b is less than 4 in most of the cases, there is one excep- 
tion: b ~ 5.6 was reported for a quasi-two-dimensional 
(2D) organic material of the BEDT-TTF family^ A de- 
viation of b from 4 can also be inferred from the optical 
data on UPt3ji£ Sr 2 Ru04r^and Crji£ see Ref.[8|for more 
details. 

Motivated by these findings, we revisit the fi/T scal- 
ing of the optical conductivity of a FL in this paper. 
We extend the Eliashberg's analysis of the Kubo formula 
for the conductivity to finite Jl and obtain an expression 
for er(f2, T) to all orders in the electron-electron inter- 
action. To discuss the results, it is convenient to iden- 
tify two distinct frequency regimes, the "high-frequency" 
and "low-frequency" ones, and also two types of FLs, the 
"conventional" and "non-conventional" ones. 

As far as the frequency regimes are concerned, f2 is 
larger than IraE R (Q,T) in the high-frequency regime, 
while £1 < lmY. R (il, T) in the low- frequency one. (The 
low- frequency regime also includes the dc limit of fl = 0.) 
In the context of the standard Drude formula, these 
regimes are also referred to as "non-dissipative" or "re- 
active" and "dissipative" , correspondingly. 

Turning to two types of FLs, we define a "conven- 
tional FL" as such in which the leading to and T de- 
pendencies of lm'S R (uj,T) are given by lj 2 + ir 2 T 2 , as 
in Eq. while the higher-order terms may be non- 

analytic. In a "nonconventional FL" , already the leading 
term in lmT, R (uj,T) is a non-analytic function of ui and 
T. For a wide class of interactions that remain finite 
at q = 0, the demarcation line between the two types 
of FLs is determined by the dimensionality: the case of 
D > 2 corresponds to conventional FLs, while the case of 
1 < D < 2 corresponds to non-conventional FLs. In the 
latter case, ImS H (a;,0) oc \u\ D and ImS H (0,T) oc T D . 
In the marginal case of D = 2, ImE R (uj,0) cx w 2 ln|w| 
and ImS fl (0,T) cx T 2 lnT. 

In the high-frequency regime, current-carrying quasi- 
particles can be considered as nearly free, so that the 
residual interaction among quasiparticlcs, which gives 
rise to their finite lifetime, can be treated as a pertur- 
bation. We show that in this regime Rep(Jl, T) is given 
by Eq. (|1.2p for both conventional and non-conventional 
FLs, as well as for the marginal case of D = 2, despite 
qualitative differences in the self-energies in these cases. 
Our analysis keeps full track of the vertex corrections to 
the conductivity and thus takes both normal and Umk- 
lapp scattering processes into account. We argue that the 
Air 2 coefficient of the T 2 term in this formula is a conse- 
quence of the "bosonic first-Matsubara- frequency rule", 
which stipulates that a function obtained by analytic 
continuation of Rep(f£, T) to the first (non-zero) bosonic 
Matsubara frequency, — > ±2i7rT, does not have a T 2 
term. Equation (|1.2p obviously obeys this rule. 

In the low-frequency regime, Rep(Sl,T) differs from 
Eq. (|1.2[) because the interaction among quasiparticles 



can no longer be treated as a perturbation, and this 
affects the T 2 and SI 2 terms in Rep(il, T) in different 
ways. We analyzed the change in the functional form of 
p(Q,T) between the high- and low-frequency regimes in 
the "zero-bubble approximation"— and found that the 
change is numerically quite small, i.e., the formula 

Rep(n,T) = A' (n 2 + bir 2 T 2 ) (1.3) 

with b ~ 4 remains quite accurate down to the lowest O, 
although the exact form of Rep(Sl, T) in the entire range 
of Vi is different from that in Eq. I|1.2[) . We also analyzed 
Rep(Sl,T) in the "incoherent regime", where all energy 
scales are of the same order, i.e., SI ~ T ~ ReS fl (0, T) ~ 
lm£ fl (Sl, T), and again found a good fit by the VL 2 +bir 2 T 2 
form with b ~ 4. 

Equation (|1.3[) is to be taken with some caution, be- 
cause the zero-bubble approximation neglects the cor- 
rections to the current vertex in the polarization bub- 
ble. Physically, vertex corrections differentiate between 
normal and Umklapp scattering processes. In the high- 
frequency regime, both normal and Umklapp processes 
contribute to the resistivity provided that the Fermi sur- 
face (FS) is sufficiently anisotropi c^ 17 ' 18 (a precise defi- 
nition of "sufficiently anisotropic" is given in Sec.lIV[). As 
a result, Eq. (jl.2p remains valid when the vertex correc- 
tions are included. The only change is that the prefactor 
A' now contains a sum of normal and Umklapp scattering 
amplitudes. In the low-frequency regime and, in partic- 
ular, at fi = 0, the resistivity of an impurity- free system 
is non-zero only in the presence of Umklapp scattering, 
although normal processes also contribute once Umklapp 
processes are allowed^ As a result, the prefactor A in 
the dc resistivity 

p(0,T) = in 2 AT 2 (1.4) 

contains some function of the normal and Umklapp scat- 
tering amplitudes rather than just their sum, and is 
therefore different from A' in the high-frequency limit. 
What remains to be seen is how the functional form 
of Rep(Sl,T) evolves between the dc and high-frequency 
limits beyond the zero-bubble approximation. 

We then discuss the experiment, focusing mostly on 
recent optical measurements on URu2Si2> 8 » Given that 
the observed values of the coefficient b are substan- 
tially different from the FL value b = 4, we argue that 
the existing optical data cannot be explained only by 
the electron-electron interaction. Following an analogy 
with the Kondo effect^ we propose a phenomenological 
model which, in addition to electron-electron scattering, 
contains also elastic scattering by some localized decrees 
of freedom, e.g., magnetic moments or resonant levels. 
In this model, the self-energy is a sum of two parts: the 
elastic one, described by an uj 2 term , and the inelastic 
one, described by the standard FL term, lu 2 + ir 2 T 2 , i.e., 

ImS fl (w, T) — C [aoj 2 + {lu 2 + tt 2 T 2 )] , (1.5) 

where the relative weight of the elastic and inelastic con- 
tributions, a, is an adjustable parameter of the model. 
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Within the zero-bubble approximation, the coefficient b 
in Eq. (|1.3[) is related to a via 

("I 

a + 1 

For practical purposes, the model is meaningful only for 
— 1 < a < oo; consequently, 1 < b < oo. The FL value 
of 6 = 4 is reproduced for a = 0. The opposite limit of 
a = oo (and thus 6=1) corresponds to a purely elastic 
scattering mechanism. The range 1 < b < 4 corresponds 
to a mixture of elastic and inelastic mechanisms with a > 

0, whereas b > 4 corresponds to an elastic contribution 
with -1 < a < 0. 

In this classification scheme, the value of b w 1, re- 
ported in Refs.lUIIl indicates a purely elastic scattering 
mechanism, whereas b « 2.3 (and thus a w 1 .3) and 
b R3 5.6 (and thus a ~ —0.35), reported in Refs. [ll| and 
H3 , correspondingly, point at a mixture of elastic and 
inelastic mechanisms with opposite signs of the elastic 
contribution. 

We discuss one possible mechanism that leads to b m 1 , 

1. e., scattering at resonant levels, and show that this 
mechanism explains the data on URu2Si2 reasonably 
well. We refrain from identifying the microscopic origin 
of the resonant levels (except for noting that extrinsic res- 
onant impurities can hardly be the culprits) but merely 
surmise that intrinsic deep electron states can play a role 
of incoherent resonant scatterers at relatively high ener- 
gies, where a coherent Bloch state is not formed yet. 

Whereas the resonant-level model explains the optical 
data in the "6 = 1" materials, the T dependence of the 
dc resistivity can be explained only by invoking a suffi- 
ciently strong electron-electron interaction which, when 
combined the resonant elastic scattering, does not signif- 
icantly affect the optical scattering rate. We show that 
dc and optical measurements probe different scattering 
mechanisms: while a dc measurement is sensitive to both 
elastic and inelastic mechanisms, an optical measurement 
probes primarily the elastic channel, 

The rest of the paper is organized as follows. In SecHIl 
we discuss the optical conductivity of a FL in the high- 
frequency regime within the Kubo formalism. In Sec. lIIIl 
we extend the analysis to both the low-frequency and in- 
coherent regimes within the zero-bubble approximation 
for the current-current correlator. In Sec. IIV1 we discuss 
the interplay between normal and Umklapp contributions 
to the resistivity in different frequency regimes. Section 
IVl addresses comparison to the experiment. In Sec. IV Al 
we discuss the status of the experiment and conclude that 
it cannot be explained within the model which includes 
only the electron-electron interaction. In Sec. lVBl we in- 
troduce a phenomenological model which combines elas- 
tic and inelastic scattering mechanisms, and classify the 
observed values of the coefficient b within this model. In 
Sec. IV C[ we apply the resonant-scattering model to the 
data on URu2Si2- Section IVTl presents our conclusions. 



II. OPTICAL CONDUCTIVITY 



OF A FERMI LIQUID 
A. Kubo formula: rigorous treatment 

1. Preliminaries 

As in I, we consider an electron system on a lattice. 
We assume that the FS does not have nested parts and 
is away from van Hove singularities but is otherwise ar- 
bitrary. Near the FS, the bare electron dispersion, 
(measured from the Fermi energy) is approximated by 
£ k = v k F ' s — kp), where kp is a vector in the di- 
rection of k and residing on the FS. Following the con- 
ventional FL methodology, we divide electron states into 
"low-energy" (near the Fermi energy) and "high-energy" 
ones. Effects of the interaction via high-energy states are 
parameterized by the self-energy E^(cj). An expansion 
of ££(w) near the FS 



1 



Z 



-1 +(u kF -v° F )-(k-k F ) (2.1) 



(2.2) 



defines the quasiparticle renormalization factor 

\ -l 



1 



duj 



j=0 



and renormalized dispersion e k = v kF • (k — kp) of the 
low-energy states, where 

v kF = Z kF u kF = Z kF (v kF - V k S k (0) | k=kF ) . (2.3) 



(As in I, we define the single-particle self-energy as 
G^ 1 (o;,r) = u) + X k (u;,T) — e k .) The renormalization 
factor and both velocities (v kF and u kF ) are defined at 
point kp oi the FS and, in general, vary over the FS. The 
(Matsubara) Green's function describing the low-energy 
electron states is given by 



G k K0 = 



iu} m /Z kF - u kF • (k - k F ) 



(2.4) 



where u> m = irT(2m + 1). 

The combination of properties formulated above de- 
fines the "bare" low-energy theory described by the ac- 
tion 

s = T ^2 Vw,k [iuJm/Zkp - u kF • (k - k F )] ip Um ,u, 

(2.5) 

where f, is a shorthand notation for J d D k/(2ir) D . The 
residual interaction between low-energy quasiparticles is 
described by an instantaneous potential C/ q , which is al- 
ready dressed up by high-energy states and assumed to 
be non-singular for any q that connects two points on 
the FS, including q = 0. Dynamic screening of the in- 
teraction by low-energy states, which gives rise to finite 
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FIG. 1. Diagrams for the current-current correlation function: 
JCi (diagram 1) and /C2 (diagram 2). The "four-momenta" 
in diagram 2) are K — (aj m ,k), K' — (oj m /,k'), and Qo = 
(fin, 0). The shaded box is the vertex r kk ' (tj m , u) m i , Q,n). 

lifetime of quasiparticles and hence finite conductivity, is 
treated explicitly. To avoid double-counting, we assume 
that mass renormalization is already absorbed into the 
parameters of the bare theory. 

In the presence of an external electromagnetic field de- 
scribed by vector potential A, the momentum k in the 
bare action, Eq. (|2.5[) , is replaced by k — eA/c. The cor- 
responding current vertex contains the "charge velocity" 
u kF [Eq. (|2.3[) ]. Note that u kF is renormalized only by 
the k dependent part of £ k (cj);£i2i in contrast to the full 
Fermi velocity v kF , which is renormalized both by V k £^ 
and dt R p /duj. 



K 2 (n n ,T) = -^-T 2 J2 ! ! u kF • u k , F G k ( Wm )G, 



To simplify notations, we assume that a metal has cu- 
bic symmetry, in which case the conductivity tensor re- 
duces to <xy = Sijcr. The diagonal component of the 
conductivity is given by the Kubo formula 

HQ, T) = ^ [)C R (n, T) + /C dia (T)] , (2.6) 

where K, R {Sl,T) = fcf (Q,T) + K, R {VL,T) is the retarded 
current-current correlation function, represented by a 
sum of two diagrams in Fig. [TJ and (e 2 /ifl)IC dla (T) is 
the diamagnetic part of the conductivity, which cancels 
the n = term in IC R (n, T) [the sum JC R {{), T)+/C dm (T) 
must vanish for a normal metal by gauge invariance] . In 
what follows, we assume that the O = piece is already 
subtracted from K. R (Q, T) and do not specify an explicit 
form of /C dla . 

On the Matsubara axis, the diagrams in Fig. Q] are 
given by 

/Ci(0„,T) = -jjT^2 [ i& F G k (u m )G k (u; m + O n ) 
w m " /k 

(2.7) 

and 

(uj m + Sl„)r kk '(a;„ l , w m ', fi n )G k '(a; m ')G k '(a; m ' + Q n ),(2.8) 



where the Green's functions and the vertex part 
r kk '(wm, kW, Q n ) contain the effects of residual inter- 
action between low-energy quasiparticles. 

The calculation of diagram 1 in Fig. [T]is fairly straight- 
forward. Replacing the Matsubara sum by a contour in- 
tegral and converting the momentum integral into inte- 
grals over de k and over the FS element eL4 kF , we obtain 
for the imaginary part of K, R 

I m /Cf(Q , T) = fdA^ J / cfc, / d£k 

x [n F (uj) - n F (uj + ti\ ImG R (w)lmG R (u) + O), (2.9) 

where np(e) is the Fermi function, G R ' A (oj) = Z-^ F /[ui — 
e k ±iZ kF ImS k (w)], and ImE k (lu) accounts for the resid- 
ual interaction. 



given by Eq. (jl.ip with the only proviso that the prefactor 
now varies along the FS: G — > G kF . The dependence of 
ImE kfc (tj,T) on e k is weak and can be neglected. The 
integral over e k is then solved readily: 

J de k lmG R (Lu)lmG R (Lu + SI) (2.10) 
= 7rZ kF Im [n/Z kF + zImE« F («, T) + ilmS^ (w + SI, T)] ~* . 

The high-frequency regime is defined by the condition 
ft»Z kF Im£ kF (r!,T). (2.11) 



2. Canonical Fermi liquids 

In this Section, we restrict the analysis to conventional 
FLs. (We will show later, in Sec. Ill A 3l that the result for 
the conductivity applies to non-conventional FLs as well) . 
For a conventional FL on the lattice, ImS ki7 (lu, T) is still 



For a conventional FL, this condition implies that SI 3> 
Gk F max{ri 2 ,T 2 } for all points on the FS. The relation 
between fi and T is arbitrary but we do assume that 
Sl,T <C Ef- In this regime, Eq. (|2.10[) is expanded in 
the imaginary parts of the self-energies and their sum is 
averaged with the difference of the Fermi functions. For 
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a conventional FL, the last step amounts to 

du [n F (uj) -n F (u) + Sl)] [imE^ (w) + ImEjj^ (u> + Q)] 



= C* kF / <£o;[ra F (cj)-ra F (cj + ft)] [w 2 + (w + O) 2 + 2(ttT) 2 



(2.12 



-c kF n (n 2 + 4tt 2 t 2 ) 



It is at this step when the difference between the coef- 
ficients of the T 2 parts in ImE^ and a occurs. Using 
flUEl) , we obtain 



Reeri(ft,T) = — Im/Cf (fi,T) = Si 



r> 2 + 4^ 2 T 2 
S12 



(2.13) 

with B x = (4e 2 /3B(2^) D )^dA kF ( U 2 F /« kF )Z3 F C kF . 
In the high-frequency regime, Reer <C Imcr = uj 2 /AitVL, 
where uj p is the effective plasma frequency. Expanding 
p(0,,T) = l/a(Q,T) in Recr/Imer, we obtain Eq. (TO)) 
with prefactor A' = (4tt) 2 Bi/uj*. 



To analyze the contribution of the vertex corrections 
represented by diagram 2 in Fig. [TJ we perform ana- 
lytic continuation of K2iSln.1T), following the procedure 
developed by Eliashberg- 5 The resulting expression is 
quite involved but to find the real part the conductiv- 
ity we need only that part of K,^(il,T) which contains 
the product of the retarded and advanced Green's func- 
tions located on the same side relative to the vertex. 
Only such products will survive upon integrating over 
e k and e k '. In general, K,tf(£l,T) contains vertices which 
are obtained by analytically continuation the Matsub- 
ara vertex r kk < (w m , u) m i , Q n ) via the following relations: 
iio m = w+ilmw, iu m i = w' + ilmw', and iVL n = tt + ilmQ, 
where all the imaginary parts are infinitesimally small. 
Analytic properties of continued vertices are determined 
by relations between the imaginary parts of the three fre- 
quencies. The part of Kf{£l, T) that we are interested in 
contains vertices rJ i I k ^ IV (oj, uj', tt), where Roman numer- 
als indicate regions in the (Imw, Ima/) plane, as shown 
in Fig. [3] (for definiteness, we set ImSl > 0). Explicitly, 



/Cf(fi,T)= j J u kF u K JdujJ dJ Mu)-n f (w + n)] Gg(u + n)G£(uj)T Uk ,(u,u',n)G$iu' + n)G&(uj'), 

(2.14) 

where 



r^^coth^^-r™ 



coth 



LO+Uj' + Q 

2T 



0, 



rS' - r k V k - tanh ^r", + tanh ^r k v k , (2.15) 



2T 




K K' 

FIG. 2. Lowest order diagrams for the vertex in diagram 2) 
of Fig. ED 



[For brevity, we do not spell out the arguments uj, uj', O 
which are the same in all vertices in Eq. (j2,15[) ,] 



Equations (|2. 15|) and ()2.14p allow one to extract the f2 
and T dependences for any vertex diagram. For example, 
vertex diagram a in Fig. [5] reads 

r^(w m ,w TO ',n n ) = t/k-k' n k-k'(w m - w m /),(2.i6) 

where II q (e) is the polarization bubble. Continuing 
this expression to real frequencies, we obtain T^, 1 = 
t^ k _ k /IIk_ k / (u) — uj') in region II, where Im(w — uj') > 

o, and r» k n = = uZ_ v n£_ v ( u ~ w ') = 

^k-k'Q 2 [n k _ k /(w — uj')] in regions III and IV, where 
Im(w — uj') < 0. Combining the contributions from re- 
gions II-IV, we obtain 



r kk' = V£-v { 2ReII k-k'(w ~ w ') [»f(w') ~ «f(w' + 0)] + 2iImIT k i _ k ,(w - uj') [2n B (u' - w) + n F (a/) + n F (w' + Q)]} , 

(2.17) 

I 



where ris(e) is a Bose function. 

As before, we replace each of the two momentum in- 
tegrals in Eq. (|2.14p by integrals over the Fermi surface 



and over the dispersion, and set k = k F and k' = k' F 
everywhere except for the Green's functions In the high- 
frequency regime, the Green's functions in Eq. (|2.14l) can 
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Imco 




Im (o>K0H-£2)=O 
FIG. 3. (color on-line). Regions of the (Ima;, Imu/) plane. 



be replaced by the bare ones [Eq. (|2.4I) ]; then the prod- 
uct fde k Gg{uj + Sl)G£(u) f devGSiu 1 + fi)Gj*(w') = 
—Aw 2 Z'^ f Z^, /ft 2 is real. Therefore, the imaginary part 

of the current-current correlator is given by ImF kk l from 



Eq. (pOTl) . Recalling that ImII q (0) 



kk' 

-T>qfl, where 



(27T) 



dA 



k' 



Zk' Zk' +q*( e k' +q) e k ,=0 



(2.18) 



(cf. Eq. (2.7) of I), and relabeling ui' — » us + fl', we obtain 
the contribution of diagram a to the conductivity 



{a} 



Reer. 



where 



(2tt)* 



gL4l 



k' 



da; 



The integral over fi' in Eq. (|2 . 19[) is the same as in the 
sum of the imaginary parts of the self-energies corre- 
sponding to diagram a in Fig. 1 of I [cf. Eq. (2.5a) in 
I]. This integral gives \ {uj 2 + (u> + ft) 2 + 2n 2 T 2 ). Av- 
eraging the last result with the difference of the Fermi 
functions, as in Eq. (|2.12p . we obtain 



u k' 



dn'Af{uj, n', o) 

-k', (2-19) 



x [2n B (n' 



n F [oS) — n F (uj + f2)] 



doj / dfl'Af{u,n',Q) = 



i n 2 + 47r 2 r 2 

3 ft 2 



(2.21) 



Thus Recr^ differs from Reci in Eq. (|2. 13|) only by a 
prcfactor, which can be read off from Eq. (|2.19j) . 

Other diagrams for JC^(Q,,T) can be analyzed in a 
similar fashion. For example, diagrams b and c in 
Fig. [5] are similar to diagram a with the only difference 
that Ul_ k ,Imng_ k ,(uj - uj') in Eq. (|2TT7| is replaced by 
^k-k'InxPj^, (uj' — uj), where 



Im7>£ k ,(fi) = f [n F (e)-n F (e 



n)] 



(2.20) 

n)] . to 
j 



x / ImG p (e)ImGp +k '_k(e + 0)C/k- P . (2.22) 

P B 

In I, we showed that the analytic properties of V in the 
frequency plane are the same as those of the polarization 
bubble. Therefore, the contributions of diagrams b and 
c to Rec, which are equal to each other, also scale as 
(Q 2 + 47r 2 T 2 )/f2 2 with a prefactor different from that of 
diagram o. [For J7 q = const, the combined contribution 
of b and c cancels that of a.] 

The Cooper-channel vertex-diagram d-appears to be 
somewhat different from particle-hole diagrams a-c but, 
in fact, it gives the same result. To see this, we notice 
that the Matsubara vertex in diagram d depends only on 
the combination uj m + u) m i + fl n ; hence, analytic prop- 
erties of the retarded vertex depend on whether one is 
above or below the Imuj + Imu/ + ImO = line in Fig. [3] 



Therefore, r£ k , = r kk , = (rj^,)", and the vertex reduces 



r kk' = 2ReC k \,(w, uj', Q) [n F (uj') - n F {uj' + ft)] + 2TmC k 7 k ,(w, uj' , 0) [2n B (u' + w + Q) + n F {uj') + n F (uj' + O)] 



(2.23) 



As before, we need only the imaginary part of of the 
vertex which contains 

ImC&, =JI ^ImG*(e)lmG* +k ,_ p (uj W + O - e) 



x tanh — (7 k _pC/ p _ k ' . 



(2.24) 



[For J7 k = const, C reduces to a Cooper bubble.] Sub- 
stituting explicit expressions for the spectral functions 
and integrating over e and e p , we obtain ImC^, — 
Cq(uj + uj' + Q) in the low- frequency limit, where Co con- 
tains a product of two interactions averaged over the 



I 

FS. Substituting this result into the imaginary part of 
Eq. (|2.23p . and relabeling u> — > —uj and uj + uj' + SI — > Of, 
we again arrive at the same integral as in Eq. (|2.21j) . 

The recipe for extracting the fl 2 + 4ir 2 T 2 scaling form 
of Rep from a diagram of arbitrary order is now clear: one 
needs to extract a factor of from either a particle-hole 
or particle-particle convolutions of the Green's functions 
in the vertex, integrate it with the combination of the 
Fermi and Bose functions in Eq. (|2.20l) . and then aver- 
age the result with the difference of the Fermi function. 
Up to a prefactor, all diagrams produce the same fl/T 
scaling form of Rep given by Eq. (|1.2[) . As it was the case 
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with the self-energy considered in I, the overall prefactor 
cannot be expressed in a compact form. 

3. N on- canonical Fermi liquids 

The analysis of the preceding section was limited to 
the case of a conventional FL, when the infrared sin- 
gularity arising from the f2/g scaling of the polariza- 
tion bubble is suppressed by the phase space volume, 
and each of the diagram considered above is convergent 
on its own. In D < 2, the phase space is too small 
to suppress the singularities and the self-energy scales 
with a; in a non-canonical way: as oj 2 ln|w| in D = 2 
and as \oj\ d in 1 < D < 2. However, infrared singu- 
larities in different diagrams for the conductivity can- 
cel each other. This cancelation manifests the gauge- 
invariance of the conductivity. In perhaps more famil- 
iar terms, this effect makes the conductivity to depend 
on the transport rather than single-particle relaxation 
time. It is more convenient to see this effect in the Boltz- 
mann equation, where the electron-electron contribution 
to the conductivity is expressed via a change in the elec- 
tron current carried by two electrons before and after 
a collision: 17 - 18 (Av) 2 = ((v k + v p - v k _ q - v p+q ) 2 ), 
where (...) stands for averaging over the FS. For q — > 0, 
(Av) 2 vanishes as q 2 , which suppresses the infrared sin- 
gularity for D > 1. 

To see how the same cancelation occurs in the Kubo 
formula, we consider two diagrams: diagram 1 in Fig. [T] 
with both Green's functions dressed by a single-bubble 
self-energy correction (diagram a in Fig. 1 of I) and 
diagram 2 in Fig. Q] with vertex correction a in Fig. [3] 
As we are interested in the q = limit, it is convenient 
to decompose the momentum transfer q into components 
along and perpendicular to the local Fermi velocity: q = 



?||Vk + q±, where v k = Vk/«k, q\\ « ?i « k F , and Uf 
is the characteristic "radius" of the FS. Accordingly, the 
(renormalized) dispersion is expanded as 

£k+q = £k + «k<?|| +q±/2m k , (2.25) 

where m k measures the local curvature of the FS. The 
imaginary part of the self-energy insertions into diagram 
lof Fig.[T]) is given by [cf. I, Eq. (2.5a)] 

Im££'>,T) = Jul J ^[n B (n)+n F (u> + Q)} 

xlmG k i +q (w + f^ImH^ (O), (2.26) 

where, as before, Imn q (f2) = — 2? q £l and 2? q cx 1/q at 
<7 — !• 0. We neglect q\\ everywhere but in ImG^ +q (cj + ft), 
integrate over q\ | , and substitute the result into Eq. (|2 .9[) , 
which is then expanded in ImS fl . This yields 

Im/Cf(n,T) = A| ?0!L<Lzl p j du J <KM(w, O', 0) 

X / l^fe^V^q^+qx, (2-27) 
where Af is given by Eq. (|2.20l) . In D < 2, the qj_ integral 
is Eq. (| 2. 2 71 infrared divergent. However, this divergence 
is canceled by vertex part a in Fig. [2] To see this cancela- 
tion, we need to assume that not only the charge velocity, 
defined by Eq. (|2.3|) . but also its derivative on the FS is 
known. Then, relabeling k' = k + q in Eq. (|2.14p . we 
expand u k + q as 

u k +q = u k + ([v k <7|| + qx] • V) u k . (2.28) 

Substituting this expansion into Eq. (|2.14j) with V — 
, we obtain 



Im/C* (fi, T) = ^j (2 ^ )g k ; kF J ds k J ^ J jduj dO'uk, • (u kF + ( [v k<Z|| + qj ■ V) u k | k=kF ) 

xG k ff (w + Q)G^{u)G^(u + + f!)G k+q (w + tf)U%D*fif(u, O', O). (2.29) 



Hg G k+q u k+q over q\\ 



First, we integrate the product G k , 
setting g|| = everywhere else in the integrand. The qu 
independent and linear- in-gj | terms in Uk+ q produce two 
integrals 



2tt 



G^q (w' + fi)^ (o/) = 



VkQ 
^k F +g ± q ± 



The term proportional to qj_ vanishes by symmetry, and 
we obtain 



2m k 

(2.30) 



where terms of order f2 



T were neglected compared to 

A 



q^l/2m] l in the second line. Next, we integrate G^G k 
over e k , setting e k = everywhere else in the integrand. 



ImJC£(fi) = -^ 
u kF 



dA v 



(27T) 
w k F , 



du 



where 



iD-l, 



(2*) 



D-l 



2mk F Wkj 



7 2 
Z k F - 



(2.31) 



(vk • V k ) u k |k=k 

(2.32) 
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The sum of two contributions to Im/C, i.e., Eqs. (|2.27[) 
and (|2.3ip . contains a combination 



„ kF Z^ p Z kF+ q u 2 . F (Z ]iF - Z kF+cu ) (2.33) 



2nik F fk, 



"Uk F • (Vk • V k ) Uk |k=k F ^k F +qj 



which vanishes as q\ for q± — > and thus suppresses 
the l/q± divergence for all D > 0. Therefore, the optical 
resistivity in the high-frequency regime has the same Q 2 + 
4ir 2 T 2 form both in conventional and non-conventional 
FLs. Notice that the other transport coefficients behave 
differently in these two cases; for example, the dc thermal 
conductivity of a 2D FL behaves as 1/TlnT, as opposed 
to the l/T behavior in 3D^i 



B. Bosonic flrst-Matsubara- frequency rule 

Just as it was the case for a single-particle self-energy 
considered in I, the ft, 2 +4ir 2 T 2 scaling form of the optical 
resistivity can be related to the analytic properties of the 
current-current correlator along the Matsubara axis. 

First, we consider the Matsubara version of dia- 
gram 1 in Fig. [TJ given by Eq. (|2.9j) . Recalling that 
sgnSk F (w m ) = sgnw m , we integrate over £k to obtain 
(for Q„ > 0) 



2iT 



D(2tt 



,D-1 



cL4 k 



Uk, 



(2.34) 



E 



Wm =-n n -HrT Z kj 

For n„ = 2itT, only one term with u) m — —nT survives 
in the sum 



Kx{2ttT,T) = - 



2iT 



D{2tt) d - 1 
I 



^ + Sk F (^r) + Ek F ("7rr)- 



(2.35) 



Since the single-particle self-energy satisfies the first 
(fcrmionic) Matsubara frequency rule, i.e., £k F (7rr) = 
^k F (— kT) = + 0(T D ), the residual interaction drops 
out from JC 1 {2nT,T) (to order T 3 ). Consequently, 
Im/Cf (fi,T) vanishes (again, up to 0(T D ) terms], when 
continued to the first non-zero bosonic Matsubara- 
frequency 2iriT. 

Next, we integrate over £k and e^' as over indepen- 
dent variables in the Matsubara version of the vertex 
part, Eq. (12.81) . setting £k = £k' = in the rest of the in- 
tegrand. At n = 27rT, only the w m — uj rn i — —ttT terms 
survive in the fermionic Matsubara sums and, as before, 
the fermionic self- energies, evaluated at ±7rT, drop out. 



Therefore, 
JC 2 (2ttT,T) 



D(2tt 



dA v 



dA 



k' 



Uk F • u k ^ 

Wk F U k ' 



xZ 2 F Z 2 r kFjkJ? (-ttT, -ttT, 2ttT). 



(2.36) 



Particle-hole diagrams for ^^(WmjWm/, O n ), e.g., dia- 
grams a-c in Fig. [2j depend on uj m — Lo m > , while particle- 
particle diagrams, e.g., diagram d, depend on u) m +oj' m + 
£l n . Since both combinations of the frequencies vanish 
at uj m — u) m i = —ttT and fl n = 2ttT, the vertex in 
Eq. (|2.36p is static, i.e., it does not contribute to the real 
part of the conductivity. Therefore, the vertex-part con- 
tribution to lmK, R (il, T) vanishes as ft = 2mT as well. 

Strictly speaking, the proof presented above is 
valid only for conventional FLs, because the single- 
particle self-energy obeys the fermionic first-Matsubara- 
frequency rule only in this case. However, as it was the 
case in Sec. Iff A~3l deviations from the canonical behavior 
caused by infrared singularities must cancel between dif- 
ferent diagrams. We did not attempt to repeat the proof 
for non-conventional FLs because the final result, ap- 
plicable to both conventional and non-conventional FLs, 
clearly shows that Reer vanishes at Vl = 2iirT. 



III. KUBO FORMULA: 
ZERO-BUBBLE APPROXIMATION 

A. High-frequency regime 

Having shown that all diagrams for the conductiv- 
ity produce the same scaling form in the high-frequency 
regime, we now consider the case of lower frequencies, 
when Q < Z kF S kF (^,T). The full analysis of the Kubo 
formula in this regime is rather involved. We will simplify 
our task and focus on diagram f in Fig. [T] which does not 
include vertex corrections. Although such an approxima- 
tion can be rigorously justified only in a few special cases, 
e.g., in the D — oo limit of the Hubbard model^ or when 
electron-electron scattering connects points of FS with 
mutually perpendicular Fermi velocities^ it provides a 
convenient way to describe a crossover between high- and 
low- frequency regimes. We also adopt a slightly different 
version of the FL theory, compared to that considered 
in the preceding part of the paper. Namely, we assume 
that the self-energy is isotropic and local, i.e., that it 
depends on w much stronger than on k — Icf, and in- 
clude the effects of interactions at all energy scales into 
the self-energy which, to order w 2 ,T 2 is now given by 
E fl (w, T) = ui(l + A) + iC (uj 2 + n 2 T 2 ) . (Since we ne- 
glected the variation of the self-energy over the Fermi 
surface, the subscript will be suppressed from now 
on.) Since this model accounts for effects on interaction 
at all energy scales, the charge and Fermi velocities en- 
tering the conductivity diagram coincide with the bare 
Fermi velocity: u kF = v kF = v£ f . Integrating over the 
bare dispersion, e k , we obtain for the conductivity given 
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by diagram 1 



B. Low- frequency regime 



u p0 



4ttS1 



dco- 



rijr(u — f2) — uf(u;) 



+ £ fl (w, T) + T, R {n — u),T)' 
where the bare plasma frequency is given by 



/<<> 



2e 2 



4tt L>(2tt) 



(3.2) 



To cast the high-frequency limit of Eq. (|3 . 1 1) into a form 
of the "extended Drude formula",— we expand in ImE^, 
evaluate the frequency integral, and bring the result of 
integration back into the denominator, which yields 



4tt n(l + A) + z 2 ^ (S7 2 



4tt 2 T 2 ) 



(3.3) 



Rep?^) = (SI 2 + 4. 2 T 2 ) 



LUp Q 3 



(3.4) 



where HF stands for "high frequency". We remind that 
Eqs. (|3.3|) and (I3.4[) are still valid only the high-frequency 
limit, defined by Eq. (|2TTT]) . 



We now analyze (7i(S7, T) at SI — > and finite T, when 
condition (12. lip is no longer valid. 

In what follows, we will need the numerical values of 
integrals 



dx 



1 



o cosh x (x 2 + ir 2 /A) 



n+l 



(3.5) 



in particular, I = 0.333..., I\ = 0.117..., 1% = 
0.043 . . . Substituting 51 = into Eq. (|3.ip and inegrating 
over w, we obtain the dc resistivity as 



PT F (0,T) 



47T 



(a C) 4tt 2 T 2 



u p0 



where LF stands for "low frequency" and 
2 



= 0.608. 



(3.6) 



(3.7) 



On the other hand, extrapolation of the high-frequency 
conductivity in Eq. (|3.3|) to SI = gives 



(O4 0.T) 



4tt 2C 

Wp 3 



4tt 2 T 2 



(3.8) 



We see that, in the zero-bubble approximation, the pref- 
actors in Eqs. (|3.4|) and (|3.8j) turn out to be very close to 
each other: (2/3)/0.608 = 1.097. 

We can also obtain the frequency dependence of 
<7i(f2, T) at SI — > by expanding Eq. (|3.1|) further in 
SI and casting the result into the form of Eq. (|3.3[) . Ex- 
panding in Q and evaluating the integrals over w, we 
obtain 



^ F (Q,T)=i 



47T 



p^ F (Sl, T) = -^-oqC ( Att'^T 2 + a 2 n z + a 3 



ifi(l + A) + ia C (4tt 2 T 2 + a 2 Sl 2 + ias ^ffi" 

(1 + A) 2 S1 2 



4tt 



u p0 



C 2 T 2 



In Eqs. p.9p and (|3.10p . ao is the same as in Eq. (|3.7p while 



01 = apm 4 '/4= 1.053. 
- 4 



«2 



«3 



a ir 



4 

32 



cosh x ( x 2 + 2^ 
Ja _ tt^oo^ =0 ()(j:;j) 



1 — 2xtanha; - 



2 1 - 2sinh 2 (a;) 

3 cosh 2 x 



r 



(3.9) 



(3.10) 



2x 2 



1.030. 



(3.11) 



In the FL regime, the imaginary part of the self-energy, ~ CT 2 , must be much smaller than T. Therefore, an ex- 
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pansion in SI should be in powers of Jl/CT 2 . This is how 
the first and the last terms (with coefficients ai and 0,3, 
correspondingly) in the denominator of Eq. (|3.9|) were 
obtained. However, because 03 happens to be numeri- 
cally very small, we also included the leading term from 
the expansion in Q/T (with coefficient 02). In practice, 
the last term in Eq. (|3.9p can be ignored so that 

47T 

Rep\ F (n, T) w —a Q C (4tt 2 T 2 + a 2 n 2 ) , (3.12) 

which is again very close to the high-frequency form, Eq. 
(EH). 

Nevertheless, a change in the ratio of the pref actors 
in the f2 2 and T 2 terms between the low- and high- 
frequency regimes indicates that the actual dependence 
of Rcpi (f2, T) is actually more complex than just a sum of 
the ft 2 and T 2 terms. We computed <7i(f2, T) numerically 
and found that the f2 and T dependences of Repi(fi,T) 
in the entire range fi,T -C Ep are well described by an 
approximate relation 



Repi( )T ) = i^ 

LU p0 6 



o 2 



3.65tt 2 T 2 ) 



(3.13) 



We see that the ratio of the 7r 2 T 2 and tt 2 terms in 
Re/?i(f2,T) is not equal to 4, but numerically is quite 
close to 4. Notice, however, that a remarkable agree- 
ment between the low- and high frequency limits is valid 
only within the zero-bubble approximation. We discuss 
effects not captured by this approximation in Sec. HVl 



C. Incoherent regime 



Equation f|3.4|) is valid in the high-frequency regime, 
as specified by Eq. (|2.11|) . Such a regime always exist in 
a coherent FL, where lmE R F (fl, T) < max{H, T}. How- 
ever, the optical conductivity of strongly correlated met- 
als is often measured in the incoherent regime where all 
energy scales are comparable, i.e., ft ~ T ~ Re£ fl ~ 
Im£ fl . Having this in mind, it is instructive to study the 
behavior of Repi(f2, T) in the incoherent regime. In gen- 
eral, calculations in this regime require a detailed knowl- 
edge of the electron-electron interaction at all energy 
scales. We use here a simple model in which ImS^fw, T) 
is assumed to follow the FL form C (uj 2 + n 2 T 2 ) all the 
way up to some cutoff frequency A and to vanish at larger 
frequencies. The KK transformation then yields 



ReS K (uj,T) = (uj 2 

7T V 



n 2 T 2 ) In 



A 



7T 



2 +7T 2 T 2 

A 2 



|A-«| 
(3.14) 



where A = 2CA/tt. At the last step, we replaced the ac- 
tual ReY, R (uj,T) by an interpolation formula which de- 
scribes the limits of both small and large (compared to 
A) frequencies but does not have a kink at uj = A. These 



forms of Re£ fl and lmT, R are substituted into the Kubo 
formula for the conductivity, Eq. (|3.1j) . and the integral 
over uj is calculated numerically. 

We found that in a wide range of and T, including 
Q ru T ~ A, the optical resistivity can be well approxi- 
mated by 



Repi(fi,T) 



4tt 



0.645 (n 2 



3.86tt 2 T 2 ) 



(3.15) 



u p0 



We see the same trend as we found earlier in the low- 
frequency regime: the f2 and T dependencies of the op- 
tical resistivity are well approximated by the tt 2 and T 2 
forms, although the actual function is more complex than 
just the sum of these two terms, and the ratio of the 7r 2 T 2 
and £1 2 terms is smaller than 4 but not far from 4. 



IV. 



UMKLAPP PROCESSES 



In SecHH we showed that any diagram for the conduc- 
tivity produces the same fl/T scaling form as indicated 
in Eq. (11.21) , with a prefactor which depends on the elec- 
tron spectrum. Since no restrictions were imposed on the 
change in the electron quasimomentum due to the inter- 
action, both the normal and Umklapp processes were im- 
plicitly taken into account. The interplay between these 
two types of processes is different, however, in different 
frequency regimes. 

In the high-frequency regime, as specified by 
Eq. (|2.11l) , the resistivity is finite already in the presence 
of only normal processes, provided that Galilean invari- 
ance is broken by a lattice. Even on a lattice, however, 
the leading, ft 2 + 4ir 2 T 2 term vanishes in several spe- 
cial cases, e.g., for a quadratic or isotropic FS in any D, 
and for a convex and simply-connected in 2Di 17 ' 18 i 23 In 
these cases, the optical resistivity scales as max{f2 4 , T 4 }. 
[The case of an isotropic and quadratic spectrum corre- 
sponds to a Galilean-invariant FL, the conductivity of 
which retains a free-electron Drude form regardless of 
the electron-electron interaction.] In what follows, we 
assume that the FS does not belong to any of the types 
specified above, so that normal processes do contribute 
to the leading term in the resistivity. If Umklapp pro- 
cesses are also allowed, they affect the resistivity as well. 
The prefactor A' in Eq. (jl.2p is proportional to the inter- 
action vertex, T. In the high-frequency regime, T is just 
a sum of the vertices for normal and Umklapp processes 
(Fn and Tu, correspondingly), i.e., 



A' ocT' 



(4.1) 



In the opposite limit of 51 = 0, the resistivity of an 
impurity-free system is non-zero only in the presence of 
Umklapp scattering. However, once Umklapp processes 
are allowed, normal processes contribute as wellj 1 ^ at 
least as a correction to the Umklapp contribution (again, 
if the FS is not of one of the types specified in the 
preceding paragraph). The effective vertex T, entering 
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Eq. (jl.4p . is now a non- trivial function of Tn and T\j 
which can be represented in the following scaling form 

A oc r = ru$(r N /ru). (4.2) 

On general grounds, one can infer that $(x — >• 0) — C\ + 
0(x) and $(a; — > oo) = C2, where Ci^ are constants. 
The ratio Tn/Tjj and the function itself depend on 
the details of both the bandstructure and the interaction 
and are by no means universal. Therefore, prefactors A' 
and A differ by a non-universal factor, which is expected 
to be of order one but not specifically close to 1. 

Even if, for some reason, normal processes are absent, 
A and A' still differ because, when calculating the opti- 
cal resistivity in the high-frequency regime, one expands 
the Green's functions in the self-energy and averages the 
result with the difference of the Fermi functions, while 
in the low-frequency regime the self-energy must be kept 
in the denominators of the Green's functions. Although 
it turns out that A and A' almost coincide in the zero- 
bubble approximation (cf. Sec. MI Al) . there is no guaran- 
tee that this remains true if vertex corrections are taken 
into account. 

We conclude this section with a remark in regard to a 
statement by Rosch and Howell,— who argued that the 
coefficients o.q and j3o in Rep(f2, T) = c(q£1 2 + /3qT 2 are 
not, in general, related. For reasons explained above, 
this statement is correct if Re/?(f2, T) is supposed to de- 
scribe the whole range of frequencies: from low to high. 
However, as we have already emphasized, the formula 
ao^ 2 + floT 2 with constant ao and /3q does not describe 
a crossover between the high- and low-frequency regimes, 
Nevertheless, a and (3 are universally related in the high- 
frequency regime, where Po/ao = Air. 

This section concludes our analysis of the conductivity 
of a FL. To summarize, we have shown that the scaling 
form of the optical resistivity in Eq. (11.2[) is quite robust. 
In the high-frequency regime, this form is produced by all 
diagrams for the conductivity. If vertex corrections are 
neglected, then one can go beyond the high-frequency 
regime. It turns out that, with only small changes in the 
numerical coefficients, Eq. (|1.2[) form works well beyond 
its nominal region of validity, i.e., both near the dc limit 
and at such high and T that the FL picture itself is not 
applicable. In other words, if the f2 and T dependencies 
of the resistivity are determined by the electron-electron 
interaction, it is impossible to avoid the FL scaling form 
with a coefficient of the T 2 term either equal or very close 
to 47T 2 . As discussed in the next Section, this is not what 
the experiment shows. 

V. COMPARISON TO EXPERIMENT 

A. Summary of experimental observations: 
Disagreement with the Fermi-liquid theory 

Now we turn to the discussion of the existing exper- 
imental data on the Cl/T scaling of the optical resis- 



tivity. Although the f2 2 dependence of Rep(f2, T) was 
convincingly demonstrated in "weakly correlated metals" 
(Au,Ag, and Cu)^ the T dependence, if measured, was 
found to result from the electron-phonon rather than the 
electron-electron interaction. This is not surprising since 
the electron-electron interaction in these metals is rela- 
tively weak and one needs to go to very low temperatures 
to observe the T 2 dependence. To the best of our knowl- 
edge, the f2 2 + 4ir 2 T 2 scaling still has not been verified 
in weakly correlated metals. 

On the other hand, the fi/T scaling of the optical con- 
ductivity in strongly correlated metals has been studied 
quite extensively; a detailed summary of experimental 
observations can be found in Ref. [|| The conclusion of 
these studies is quite surprising: fitting the measured 
"optical scattering rate" into a phenomenological form 

— = %Rep(fl,T) = const x In 2 + 6tt 2 T 2 ) (5.1) 
r(fi,T) lj£ 

has not produced b close to 4 in any of the cases studied so 
far. In some cases, e.g., in the heavy-fermion compound 
URu2Si2 8 above the 17.5 K transition into the "hidden- 
order" (HO) state and in rare-earth based doped Mott 
insulators (Ceo.095Cao.05TiO3.04 2 and Ndo.gosTiO^), b 
has been found to be close to 1 rather than to 4. On the 
other hand, a recent studyii of the underdoped cuprate 
HgBa2Cu04+a reports b ps 2.3, while another study of 
an organic material from the BEDT-TTF family reports 
6 « 5.6.1* 

In the preceding sections, we showed that b 4 is a ro- 
bust property of FLs with electron-electron interaction. 
We must then conclude that even though l/r(f2, T) in 
the compounds mentioned above exhibits FL-like depen- 
dences on and T, the lack of a FL-like O/T scaling 
indicates that these dependencies do not come only from 
the electron-electron interaction. In the remainder of this 
Section, we attempt to explain the discrepancy between 
the FL theory and the experiment. 

B. Elastic vs inelastic contributions to the 
single-particle self-energy 

In this Section, we try to identify a mechanism respon- 
sible for deviation of the observed coefficient b from the 
FL value of 4. In the preceding Sections, we analyzed the 
conductivity of a FL under an implicit assumption that 
the only scattering mechanism is the electron-electron 
interaction among itinerant electrons. However, the FL 
of itinerant electrons is not the only example of a FL. 
Another example is a FL state formed around magnetic 
impurities at energies below the Kondo temperature. In 
the Kondo case, there are two channels of interaction: an 
elastic one, which contributes an uj 2 term to the imagi- 
nary part of the self-energy, and an inelastic or electron- 
electron one, which contributes an uj 2 + ir 2 T 2 term. The 
relative weight of these two contributions depends on 
the strength of the on-site electron-electron interaction, 
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which can be conveniently parameterized by the Wilson 
ratio, R. 20 In the unitary limit, when R = 2, the elastic 
channel is twice more efficient than the inelastic one, i.e., 



lmE R (uj,T) = B--C 
o 

= B-C 



1 



n 2 T 2 ] 



-TT 2 T 2 



(5.2) 



where B is the w-independent part of the elastic contri- 
bution and C > 0. The reduction of the T 2 contribution 
to ImEfl is reflected in the optical scattering rate, which 
is obtained, as before, by substituting Eq. (15.21) into the 
Kubo formula (I3.ip (in the zero-bubble approximation) 
and integrating over uj: 



1 



r(Q,T) 



B 



2C 



(fl 2 + 2tt 2 T 2 



(5.3) 



Thus the Kondo FL belongs to a different universality 
class with b = 2. This does not explain the experiment 
yet because of the non-metallic signs of the f2 and T 
dependences of l/r(0, T) in Eq. (|5.3p . as opposed to the 
metallic signs observed in the experiment at least at the 
lowest frequencies. However, this gives us an idea to ask: 
how does a reduction of the inelastic contribution to the 
self-energy affect the relative weight of the il 2 and T 2 
terms in the optical conductivity? 

To answer this question, we introduce a phenomeno- 
logical form of the self-energy 



ImS s (w, T) = S c i(w) + C (uj 2 + tt 2 T 2 



(5.4) 



The first term describes a contribution of the elastic chan- 
nel which arises from the energy dependence of the effec- 
tive scattering cross-section. However, since scattering is 
elastic, its cross-section does not depend on the temper- 
ature (provided that the number and other properties of 
the scattering centers do not vary with T) and S e i(w) 
is T-independent. The second term describes the contri- 
bution of inelastic electron-electron interaction, which is 
the same as in a conventional FL. A particular form of 
E e i (a;) is important for determining the actual behavior 
of the optical conductivity, especially if E e i(cj) is a sharp 
function of oj, as it is the case for resonant scattering, 
considered in the next Section. For the time being, how- 
ever, we assume only that E e /(w) is an analytic function 
of uj and expand it to second order in uj as 



Eoi(w) = E c; (0) + E e ;(Q)w + aCuj 2 



(5.5) 



where the constant C [the same as in Eq. (|5.4j was fac- 
tored out for convenience, and a is another constant 
which can be of either sign. We call the elastic con- 
tribution "metallic" if a > and "non- metallic" if a < 0. 
On the other hand, the inelastic contribution is always 
metallic because C > (which is not the case for the 
Kondo model). Combining Eqs. (|5.4[> and (|5.5p . we ob- 
tain 

ImE fl -(u;,T) = E c ,(0)+E c( (0)w+C [aou 2 + (uj 2 + ir 2 T 2 )] . 

(5.6) 



The auj 2 term mimics the oj 2 dependence of the inelastic 
contribution but does not have its T 2 counterpart. We 
emphasize that the uj and T dependencies of the inelas- 
tic contribution should be consistent with the fermionic 
first-Matsubara-frequency rule, which stipulates that the 
inelastic term in Eq. (|5.4|) must vanish upon replacing uj 
by ±inT. This rule, which is obviously satisfied with our 
choice for the inelastic part, does not allow for changes in 
the relative weight of the uj 2 and T 2 terms in this part. 
Next, we substitute Eq. (|5.4p into Eq. (|3.ip . integrate 
over oj, upon which the linear-in-w term in lmE, R (uj,T) 
vanishes, and obtain the optical scattering rate as 



-77^ = - + l(* + W[tf+bir 2 T 2 ] 
t(U,T) t 3 



with 



b = 



1 



(5.7) 



(5„ 



The residual term, 1/tq, contains a contribution from 
static disorder (not considered explicitly here), and we 
absorbed E e i into this term as well. 

Now we discuss constraints imposed on the parameter 
a, and thus on b. For a < — 1, the prefactor of the second 
term in Eq. (|5.7p is negative, i.e., the f2 and T dependen- 
cies of l/r(J7,T) are non-metallic. Since this does not 
correspond to any of the experiments, we discard this 
possibility. The special case of a = —1 corresponds to 
1/t(Q,T) which depends only on T but not on f2. Dis- 
carding this possibility as well, we focus on the range 
— I < a < oo, which corresponds to 1 < b < oo. The 
FL value of b = 4 is reproduced for a = 0. The opposite 
limit of a — oo (and thus b = l) corresponds to a purely 
elastic scattering mechanism. The range l < b < 4 corre- 
sponds to a mixture of elastic and inelastic mechanisms 
with a > 0, i.e., with a metallic sign of the elastic con- 
tribution, whereas b > 4 corresponds to a non-metallic 
elastic contribution with — l < a < 0, although the f2 and 
T dependences of l/r(f2, T) in this case are still metallic. 

According to this classification scheme, the value of 
b w 1, reported in Refs. for the U, Ce, and Nd- 

based compounds, indicates a purely elastic scattering 
mechanism (a = oo). The value of b w 2.3 (and thus 
a s» 1.3), reported in Ref. [Ill for the Hg-based under- 
doped cuprate, points at a mixture of elastic and inelastic 
mechanisms with comparable weights, and with a metal- 
lic sign of the elastic contribution. Finally, b s» 5.6 (and 
thus a w —0.35), reported in Ref. fl2l for the organic mate- 
rial, also corresponds to a mixture of the two mechanisms 
but with a non-metallic sign of the elastic contribution. 

The deviation from the FL behavior is the most dra- 
matic for the b — 1 case, where it appears that the 
electron-electron interaction does not play any role. How- 
ever, this conclusion would be incorrect. In the next 
Section, we discuss one example of a purely elastic scat- 
tering mechanism, i.e., scattering from resonant levels, 
and apply this model to the URu2Si2 data. We will see 
that, while the optical conductivity can be explained by 
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resonant-level scattering alone, an explanation of the T 
dependence of the dc resistivity requires invoking a suf- 
ficiently strong electron-electron interaction. 



C. Scattering from resonant levels: 
the case of URu2Si2 

In this section, we discuss the model of purely elastic 
scattering from resonant levels, located at energy uiq away 
from the Fermi energy and of width 7. The self-energy 
in this case is given by 



lm2 a (u,T) = EelM = 



0)7 



(cj - w ) 2 + 7 2 



(5.9) 



At T — 0, the corresponding optical scattering rate is 
given by 



1 



r(O,0) 



Co 

si 



arctan 



SI — LUq 



arctan ■ 



CJ 



7 7 

(5.10) 

If the resonant level coincides with the Fermi energy, 
1/t(0, 0) is purely non-metallic, i.e., it decreases as SI 
increases. If the resonant level is away from the Fermi 
energy, 1 /r(f2, 0) is a non-monotonic function of SI with a 
maximum at SI ~ loq (see Fig. 2]). The origin of the max- 
imum is clear: as SI increases from zero to uin, the rate of 
transitions from the Fermi energy to the resonant levels 
increases. When SI becomes larger than uj , the rate de- 
creases because now the energy interval from the Fermi 
energy to the resonant level constitutes only a fraction of 
the photon energy. Expanding Eq. (|5.9|) near lo — as 



haS R (u) = C07 



1 



+ 



3^ - T , ,2 



+ Y K 2 +7 2 ) K 2 +7 2 )" 



and substituting (|5.1ip into fl3.1[) . we obtain 
C o7 (3c 2 - 7 2 ) 



(5.11) 



t(Q,T) 



const- 



sr 



tt 2 T 2 ) . (5.12) 



(The linear in u) term in Eq. (|5.1ip vanishes by parity.) 
Already for a moderately narrow level, i.e., for 7 < wov3, 
the signs of both the SI and T dependences of 1/t(Q,T) 
are metallic, and 6 = 1. The behavior of l/r(fi,T) over 
a larger range of S7 and T is obtained by substituting 
Eq. (15.91) into the Kubo formula (|3.1[) and computing 
the integral over ui numerically. The results are shown in 
Fig.[SJ To compare to the experimental data on URu2Si2 
from Ref. shown in Fig. [3 we choose uiq = 12.5 meV to 
match the position of the peak in the data. All other en- 
ergies are measured relative to luq. In Fig. 7 = 0.2cl>o, 
and the temperatures are chosen to coincide with the ab- 
solute temperatures used in the experiment (18, 22, and 
25 K). Comparing Figs. [5] and we see that the model 
reproduces the characteristic features of the data, i.e., 
a non-monotonic dependence of l/r(0, T) on SI, as well 
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FIG. 4. (color on-line). Imaginary part of the fermionic self- 
energy, Eq. (15. 9|) . (left) and optical scattering rate at T — 0, 
Eq. (|5.12[) (right) for scattering at resonant impurities. 



as an (approximate) isosbestic point at S7 sa uio, where 
l/r(0, T) apparently does not depend on T. The very 
existence of the isosbestic point implies that the data 
cannot be described by Eq. (II. 2j) with a T-independent 
prefactor A'. Nevertheless, we follow the same protocol 
as used in Ref.d, i.e., we fit the SI dependence of the com- 
puted l/r(0,T) into an SI 2 function ( shown by dashed 
lines in Fig. [5j , then fit the intercept 1/t(£! — > 0,T) into 
a T 2 function (shown in Fig.[6|), and take the ratio of the 
slopes of the T 2 and SI 2 fits. This procedure gives 6 ps 0.9, 
which is within the margin of error of the experimental 
value 6 = 1 ± 0.1. 8 

The behavior of l/r(Sl,T) in the Ce and Nd com- 
pounds fRefs. l9l andfiol. correspondingly) is qualitatively 
similar to that in URu2Si2, although the ranges of SI and 
T are drastically different. In Nd .905TiO 3 , 1/t(Q,T) 
scales as SI 2 up to about 0.1 cV, followed by a maximum 
at w 0.27 eV. The SI = intercept of 1/r scales as T 2 
over a wide temperature range: from 29 to 295 In 
Ceo.095Cao.05TiO3.04, l/r(fl,T) scales as SI 2 also up to 
about 0.1 eV, followed by a tendency to saturation; but 
the maximum is not yet revealed at the highest frequency 
measured (pa 0.14 eV). The SI = intercept also scales as 
T 2 over a wide range of temperatures. These similarities 
suggest that, despite obvious differences in composition 
and energy scales in U, Ce, and Nd compounds, the op- 
tical response in all three cases is governed by the same 
mechanism. 

Elucidation of the microscopic mechanism of resonant 
levels is beyond the scope of this work, and we make just 
a brief comment in this regard. It is very unlikely that 
clean samples studied in Refs. l8T4l0l contained consider- 
able amounts of extrinsic resonant impurities. Therefore, 
resonant states must be intrinsic to these compounds. 
We surmise that /-electrons of U, Ce, and Nd atoms, 
although arranged into a sublattice, play the role of in- 
coherent resonant levels at sufficiently high energy scales 
probed in optical measurements. 
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FIG. 5. (color on-line). Optical self-energy in the resonant- 
impurity model as a function of frequency at several temper- 
atures. Absolute values of temperatures are fixed by choosing 
cjo = 12.5 meV and 7 = 0.2<jjo- Dashed lines show SI 2 fits of 
the actual dependencies. 
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FIG. 6. (color on-line). The intercept, l/r(fi 0, T), in the 
resonant impurity model as a function of (T/ujo) 2 - 
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FIG. 7. (color on-line). Experimental results for in URu2Si2 
from Ref. @. a) Optical scattering rate 1/t(Q,T). b) The 
optical resistivity at lower frequencies from the refined reflec- 
tivity. 



comes from an intrinsic mechanism. On the contrary, 
the T dependence of the dc resistivity in the resonant- 
scattering model is purely non-metallic. Indeed, it is easy 
to see that the dc conductivity, 



D. Combined effect of the electron-electron and 
resonant-level scattering mechanisms 

Although the resonant-scattering model explains the 
results of optical measurements, this model alone can- 
not explain the temperature dependence of the dc re- 
sistivity Above the superconducting transition temper- 
ature (ss 1 K) in URu2Si2, both the a- and c-axis re- 
sistivities increase with T in a quadratic manner within 
the HO phase, exhibit a kink at HO T]r , and continue 
to increase up to about 75 K, where p a goes through 
a broad maximum whereas p c starts to saturate ! 29 ' 31 ' 32 
The slopes of the increasing parts in p a . c , both below 
and above Tjf 1 , are largely independent of the resid- 
ual resistivity, 3 - which indicates that the T-dependence 



with ImT. R (uj) from Eq. (|5.9p . increases with T as T 2 ; 
therefore, p(0,T) = l/er(0,T) decreases with T. In order 
to reproduce the metallic sign of p(0,T), at least for T 
below 75 K, one needs to bring in the inelastic electron- 
electron interaction with ImS (ui, T) given by Eq. (|l.ip . 
This seems to defy the purpose of the preceding analysis, 
as we have argued that the optical data cannot be ex- 
plained by an inelastic mechanism. It turns out however, 
that a combination of elastic and inelastic mechanisms 
explains both the dc and optical data. In the "combined" 
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FIG. 8. (color on-line), dc resistivity (arbitrary units) for a 
model form of the self-energy which combines resonant-level 
and electron-electron contributions, Eq. (|5.14[) . Temperature 
is measured in units of the resonant-level energy, ojo, which 
is also chosen to coincide with the cutoff energy A. The 
resonant-level width 7 = 0.2wq- Parameter a, defined by 
Eq. (|5.15[) . measures the relative strength of the two contri- 
butions. 



model, the total self-energy is a sum of two contributions 

Col 



Im£ fl (w,T) 



(uj - oj ) 2 +7 2 



(5.14) 



+C (uj' z + ir 2 T 2 ) F 



\J up- + ir 2 T 2 
A 



This equation is the same as we introduced in Eq. ([57 
except for now the electron-electron contribution con- 
tains a smooth cutoff function F(x) 1 defined in such a 
way that F(0) — 1 and F(x) falls off faster than 1/x 2 for 
x 3> 1. The function F(x) is chosen to reproduce a slow 
decrease of the measured a-axis resistivity at higher tem- 
peratures. (Since optical experiments probe the basal- 
plane conductivity, we focus on this case.) To minimize 
the number of free parameters, we set A = ujq. The 
relative strength of two contributions to ImY, R (uj,T) in 
Eq. f|5 . 14[) is controlled by a dimensionless parameter 



04 



(5.15) 



Larger values of a correspond to a larger electron- 
electron and smaller resonant-level contribution, and vice 
versa. Using the small-o; expansion in Eq. (|5.11|) . it 
is easy to show that dp(0,T)/dT|T^o is positive, i.e., 
"metallic", already for a > 1/4. The dependence of 
p(0, T) over the entire temperature range is obtained 
by numerical integration of Eq. (|5 . 13[) with ImE^ from 
Eq. (|5.14[) . The resulting profiles of p(0,T) are shown in 
Fig. [5] for a = 0, 1, 2. As we see, the electron-electron 




FIG. 9. (color on-line). Optical scattering rate (arbitrary 
units) for a model form of the self-energy which combines 
resonant-level and electron-electron contributions, Eq. (|5.14|) . 
Frequency is measured in units of ujo. T = 22 K, 7 = Q.2uq. 
From top to bottom: a — (red), a — 1 (blue), a — 2 (green). 



contribution leads to a qualitative change in p(0,T): a 
purely non-metallic T dependence with resonant levels 
alone (a = 0) is transformed into a curve with a maxi- 
mum (a = 1,2). The a — 2 curve is already similar to 
the measured profile of p(0,T), which increases almost 
three-fold when T is varied in between T c HO (chosen as 
the lowest temperature in Fig. [5]) and the temperature 
corresponding to a maximum resistivity 29 ' 31 ' 32 On the 
contrary, the optical resistivity is largely unaffected by 
the electron-electron contribution. Figure [9] shows the 
frequency dependence of the optical scattering rate at 
fixed temperature (= 22 K) for the same values of a 
(= 0,1,2) as in the dc case (Fig. [8]). As it is obvious 
from the figure, 1/t(0,T) is practically the same for all 
three values of a, except for a small overall shift. Repeat- 
ing the same procedure as was applied to the numerical 
data in Figs. [S] and O we again arrive at the result that 
coefficient b in Eq. (|5.ip is very close to 1. 

The results presented above indicate that the electron- 
electron and resonant-level contributions to the self- 
energy affect different parts of the frequency range: 
whereas the electron-electron contribution is largely re- 
sponsible for the T dependence of the dc resistivity and 
has practically no effect on the high-frequency optical re- 
sistivity, the resonant-level contribution determines the 
optical resistivity but plays only a secondary role in con- 
trolling the dc resistivity. This happens because dc and 
optical measurements probe different parts of the elec- 
tron spectrum (cf. Fig. [T0|) . At sufficiently low temper- 
atures, i.e., &t T ^ ujq, & dc measurement probes the 
spectrum in the region u> ~ T <C ljq, where both contri- 
butions to the self-energy vary smoothly with ui (as lo 2 ). 
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FIG. 10. (color on-line). Resonant-level (a sharply peaked 
curve) and electron-electron contributions to the imaginary 
part of the self-energy. The electron-electron parts correspond 
to a = 1, 2 and were multiplied by a factor of 8 for clarity. 



The resulting T dependence of p(0, T) is just a sum of two 
T 2 terms with opposite signs, and the electron-electron 
contribution wins this competition rather easily. On the 
other hand, the optical scattering rate is controlled by 
the region u) ~ ujq , where the resonant-level contribution 
has a sharp peak and thus dominates over the electron- 
electron one, even if the latter is strong enough to control 
the dc resistivity. 

Concluding this section, we would like to emphasize 
the importance of a sharp feature in the elastic contri- 
bution to the self-energy, E e i(w). Indeed, the classifica- 
tion scheme of different behaviors of 1/t(H,T) based on 
the magnitude and sign of the coefficient a, as defined 
by Eq. (|5.5I) . would predict substantially different val- 
ues of b for the values of the parameter a used above. 
It is easy to see that, in the resonant scattering model, 
a = 3/a which, according to Eq. (|5.8|) . implies that 
b = (4a + 3)/(a + 3). This formula gives b = 1 for 
a = 0; b = 7/4 1.75 for a = 1; and b = 11/5 « 2.2 
for a = 2. Nevertheless, fitting l/r(f2, T) curves com- 
puted with a full form of S ee rather than with its Taylor 
expansion, we obtained b w 1 in all of these cases. The 
reason for this discrepancy is that the Taylor expansion 
is not applicable near a sharp peak E ee , and one needs 
to use the classification scheme based on Eqs. (|5.7[) and 
(15.81) with certain care. 



VI. SUMMARY 

The main purpose of this paper was to highlight the 
universality of the FL result for the optical resistivity, 



Eq. (|1.2p . We showed that, within the Kubo formalism 
which takes full account of vertex corrections to the con- 
ductivity, Eq. (|1.2|) holds for an arbitrary lattice and for 
any form of the electron-electron interaction, as long as 
the system remains a FL and is away from nesting and 
Van Hove singularities. In fact, the optical resistivity 
turns out to be more universal than the single-particle 
self-energy: whereas the latter is described by the conven- 
tional form given by Eq. (jl.l[) only in canonical FLs, i.e., 
in D > 2, and deviates from this form in non-canonical 
FLs i.e., in in 1 < D < 2, the former is given by Eq. (|1.2[) 
both for canonical and non-canonical FLs. We showed 
that a particular scaling form in Eq. (jl.2[) takes its roots 
in analytic properties of the optical conductivity along 
the Matsubara axis and is consistent with the bosonic 
first-Matsubara-frequency rule. 

If a system contains not only of itinerant electrons but 
also localized degrees of freedom (magnetic moments or 
resonant levels) , the functional form of the optical resis- 
tivity changes, as specified by Eq. (|5.1[) . The magnitude 
of the coefficient b in this equation depends on the in- 
terplay between inelastic (electron-electron) and elastic 
scattering mechanisms. Completely inelastic electron- 
electron scattering corresponds to b — 4; completely elas- 
tic scattering from, e.g., resonant levels, gives b = 1; in- 
termediate cases, where elastic and inelastic channels are 
mixed, correspond to 1 < b < oo. 

As far as the existing experiments are concerned, the 
value of b = 4 has never been reported. In some cases, 
including the latest detailed study of the optical con- 
ductivity in URU2S12 (Ref. Q), the coefficient b has been 
found to be close to 1, which indicates a completely elas- 
tic scattering mechanism; a recent study of the Hg-based 
underdoped cuprate reports b « 2.3; yet another study 
of the BEDT-TTF organic material reports b « 5.6. We 
considered a simple model of scattering from resonant 
levels, and showed it is capable of reproducing the major 
features of the optical resistivity in URu2Si2 above T^. 
On the other hand, the T dependence of the dc resistivity 
can only be explained in a model which combines elas- 
tic and inelastic electron-electron scattering mechanism. 
We deliberately refrained from identifying a microscopic 
nature of resonant levels, except for stating that they are 
not likely to be extrinsic resonant impurities. More likely, 
deep / states of rare earth atoms play the role of inco- 
herent resonant scatterers at rather high energy scales 
probed in optical measurements. 

If this picture is correct, it tells us something new 
about a crossover between coherent and incoherent 
regimes in heavy-fermion materials. The conventional 
scenario of this crossover is that the only energy scale 
is the Kondo temperature (Tk)- Above Tk, localized 
magnetic moments scatter electrons incoherently, as in 
diluted Kondo alloy. Below Tk, a (heavy) FL state is 
formed and localized moments do not scatter electrons 
anymore but participate in formation of a coherent Bloch 
state. The low-energy FL state is supposed to have all 
the attributes of a standard FL, in particular, the coef- 
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ficient b must be equal to 4. This scenario is probably 
correct as long as the evolution of the system is traced 
along the temperature axis. Optical measurements add 
one more dimension: frequency. In the presence of elastic 
scattering, the variations of temperature and frequency 
do not have the same physical consequences because the 
scattering cross-section depends on the electron energy, 
and thus on the frequency of light, but not on the tem- 
perature. It appears that the crossover between the in- 
coherent and coherent regimes along the frequency axis 
contains an intermediate interval, where localized states 
scatter itinerant electrons neither as Kondo spins nor as 
screened Kondo clouds but rather as resonant levels. 

Regardless of the validity of a particular model for elas- 
tic scattering, we hope that our paper will help to rec- 
ognize the importance of the interplay between f2 and T 
dependencies in the optical data. We believe that, on 
par with much studied recently Wiedemann-Franz law 
which, if satisfied, indicates not only the FL nature of 
the ground state but also complete elasticity of the un- 
derlying scattering mechanism, systematic studies of the 
coefficient b can tell us something new about the inter- 
play between elastic and inelastic channels in strongly 



correlated electron systems. 
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